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On the basis of the nonequilibrium Green's function and nonperturbative canonical transformation 
for the local electron-phonon interaction (EPI), the quantum transport through a single-molecule 
transistor(SMT) has been investigated with a particular attention paid to the joint effect of the EPI 
and SMT-lead coupling on the spectral function and conductance. In addition to the usual EPI- 
induced renormalized effects (such as the red-shift, sharpening, and phonon-sidebands of the SMT 
level), owing to improved disentagling the electron-phonon system, it has been found that, the profile 
of the spectral function of the SMT electron is sensitive to lead chemical potentials, thus can readily 
be manipulated by tuning the bias as well as the SMT-gate voltage. As a consequence, the broken 
particle-hole symmetry in this system can be clearly recognized through the phonon sidebands in 
the spectral function. These EPI effects also manifest themselves in the nonequilibrium transport 
properties of the SMT, particularly at low temperature. 

PACS numbers: 85.65.+h, 73.63.Kv, 73.23.-b, 71.38.-k 



I. INTRODUCTION 

With the remarkable advances in the ability to fabri- 
cate the nanostructures, it is possible to explore the elec- 
tronic devices based on a single molecule, i. e. the single- 
molecular transistor (SMT)^^ which might promise to 
be potential candidates for nano-elcctronics and also pro- 
vide with an controllable tool in studying the fundamen- 
tal physics in nanometer scale. 

It has been found that the transport through the 
single-molecular transistor might exhibit profound many- 
body correlated effects at low temperature, such as the 
Coulomb blockade or the Kondo effect An im- 

portant feature for the SMT is its sensitivity to the lo- 
cal molecular vibration, which will have profound im- 
pacts on the transport properties This may also be 
true in some quantum dot(QD) systems whenever there 
is a strong coupling between the dot electron and local 
phonon mode. In addition, the SMT is usually connected 
to the outside gates via the leads. Since the applied bias 
voltage to the SMT is within nanometers, a small bias 
might cause to substantial field effect, thus the transport 
through a SMT or QD is in general a nonequilibrium 
process. 

Theoretically, many efforts have been made to- 
wards the quantum transport through the SMT 
or QD with these features taken into proper ac- 
count. Theoretical approaches developed in this 
field include the kinetic equation method^ - the rate 
equation approachjii*i£ the nonequilibrium quantum 
theory^i^4i^ii& and more recently the numerical renor- 
malization group calculation^iiSiiS In doing such in- 
vestigations, several groups have already accounted for 
the electron-phonon interaction (EPI), many-body ef- 
fects, or nonequilibrium properties ^2*22i2i However, dif- 
ferent results may be obtained by different treatments 
of these effectsii2*i& Although the numerical renormal- 



ization group method can well predict the equilib- 
rium properties of the system, it cannot be applied to 
the nonequilibrium situation directly. While regard- 
ing to the present systems, a theory capable of deal- 
ing with both nonequilibrium and strong EPI situations 
should be favorite. One candidate, the Kcldysh Green's 
function, has shown to be successful for the nonequi- 
librium systemsi 22 i 23 i 24 i 25 i 26 i 27 Within this framework 
the nonequilibrium Dyson equations are solved self- 
consistently, in which the self energies due to the EPI 
can be calculated by the standard many-body diagram- 
matic technique (22*2^ or by using the self-consistent Born 
approximation^ Such an approach can meet the present 
concerns, and can be generalized to account for the ef- 
fects arising from the phonon dynamics itself^ But, in 
general, the calculations are complicate as self-consistent 
integrals are encountered. This difficulty can be circum- 
vented by directly diagonalizing the EPI terms with a 
canonical transformationf 2 ^ 2 ^ as long as the phonons are 
in the thermodynamic equilibrium and the fluctuation 
in phonons can be neglected. By an improved treat- 
ment in the decoupling approximation as presented in the 
present work, which differs from the majority of previous 
publications; 16 ! 28 the electron can be properly decoupled 
from phonon system, and yet permits to straightforward 
solutions in a compact, analytic form. 

In this paper, by combining the nonequilibrium 
Keldysh Green's function technique with the canonical 
transformation for the electron-phonon system, we in- 
tend to study the quantum transport through the SMT 
in the presence of finite bias and strong local EPI, with 
special attention paid to the effect of the improved de- 
coupling treatment for the electron-phonon system on the 
spectral function. The paper is organized as follows. In 
section II, an Anderson-Holstein model is introduced for 
the SMT in the presence of local EPI, then our main the- 
oretical framework is described, in which the commonly 
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FIG. 1: Schematic illustration for the single-molecule transis- 
tor system. 



used procedure for decoupling the electron-phonon sys- 
tem is carefully reexamined. In section III, the numerical 
results for the spectral function of the SMT electron is 
demonstrated and discussed, including both the dressing 
effects due to the EPI and the joint effects caused by its 
interplay with the SMT-leads coupling. An interesting 
phenomenon is predicted that the profile of the spectral 
function for the SMT electron can be manipulated at 
low temperature by tuning the chemical potentials in the 
leads. Consequently, the phonon sidebands in the spec- 
tral function may be quite asymmetric with respect to 
the renormalized SMT level, indicating a broken electron- 
hole symmetry. In section IV, it is shown how the effect 
of the different decoupling approximation on the spec- 
tral function manifests itself in the conductances of the 
SMT at zero temperature and finite temperatures, which 
might yield nearly the same results at high temperature, 
but definitely different behaviors at zero temperature. 
Finally, a brief conclusion is drawn. 



II. PHYSICAL MODEL AND FORMALISM 
A. Model 

As shown in Fig^ our model system is a single elec- 
tron level in a SMT or QD, which is coupled to the local 
vibration mode as well as to two metallic leads. For the 
sake of simplicity, we shall restrict ourselves exclusively 
to the effect of the EPI and chemical potentials in two 
leads, thus ignore other factors, such as the intricacies of 
the real SMT, Coulomb interaction and spin effect. The 
model Hamiltonian can be then expressed as 



H — Hi ca ds + H p h + Hd + H T , 



(1) 



where the first two terms of the Hamiltonian represent 
the noninteracting electron gas in the leads and the local 
vibration mode of the SMT, respectively, namely, 

Hi oa d s = Ckcj^Ck, H ph = u> aJa. (2) 

keL(fl) 

Here the operator (ct) creates (annihilates) a conduc- 
tion electron with momentum k and energy et m the left 



(L) or the right (R) lead, loq is the vibrational frequency 
of the molecule, and at (a) is the phonon creation (anni- 
hilation) operator. The third term describes the cou- 
pling between the SMT electron and local phonon mode 
with strength A, 



Hr 



[ eo + A (a + at)] dM, 



(3) 



where dt (d) is the corresponding creation (annihilation) 
operator of the SMT or QD electron associated with the 
energy e . The last term Ht describes the hopping of 
electron between the SMT and leads with the tunnelling 
matrix elements denoted as Vk, 



kei(-R) 



[v k 



4d 



h.c. 



(4) 



The chemical potentials in the left and right lead arc 
denoted by (J.l(r), respectively, which arc related to the 
bias, Vbias, and the average of two lead potentials, V m ich 
through V b ia S = {^l - fJ,R)/e and V mld = (ml + Mi?)/ 2e - 
By changing the external bias and gate voltages experi- 
mentally, ^l(R) an d eo can be adjusted independently. 

In the strong EPI regime, it is appropriate to elimi- 
nate the electron-phonon coupling terms in the Hamil- 
tonian by using a nonperturbative canonical transforma- 
tion, i.e., H = e s He" s with S = ^d+d (at - a). The 

transformed Hamiltonian reads H = H p h + H c i, where 
the phonon part keeps unchanged, while the electron part 
is reshaped into 



H, 



ekC^Ck + e dtd + ^ [v k ct d + h.c. . (5) 



It is clear that due to the EPI, the energy level of SMT 
is renormalized to eo = eo — g^o, where g = (X/loq) , 
and the dressed tunnelling matrix elements are trans- 
formed into Vk = VkX, where the phonon operator X = 
cxp [— (X/coq) (at — a)] arises from the canonical trans- 
formation of the particle operator e s de~ s = dX. This 
reveals that the interaction between the SMT electron 
and the local phonon mode results in an effective phonon- 
mcdiated coupling between the SMT and the lead elec- 
trons. As in dealing with the localized polaron22iSS, here 
it is reasonable to replace the operator X with its ex- 
pectation value (X), i.e., Vk = Vk cxp [— g (N p h + 1/2)], 
where N p h is the phonon population, and can be ex- 
pressed as Nph = 1/ [cxp(/3w ) — 1] with /3 = l/ksT. 
Notice that this is an important approximation made in 
the present paper, which is valid only when the hopping 
is small compared to the EPI, i.e. Vk <C A. 



B. Formalism 

By using the Kcldysh Green's function 
tcchnique2£i2L2£, the current through an interac- 
tion region coupled to two leads can be expressed 
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J = gj- J duj { [f L ( W ) r L - / fl ( W ) r fl ] (G» - G») 

+ (r L -r fi ) G<(u;)},(6) 

where is the Fermi distribution function in the left 

(right) lead, T L ^(e) = 2irp L(R) (e) \V L{R) (e)\\ p L{R) (e) 
is the density of states in the left (right) lead, V^ R ){e) 
equals Vkei(R) f° r € = 6 k, and G r ^ a \uj), G < {uj) are the 
Fourier transformations of the standard Kcldysh retarded 
(advanced) and lesser Green functions for the dot elec- 
tron, respectively^. To note that the parameters ^l(r) 
are not the renormalizcd ones, because derivation of this 
current formula only relics on the hopping terms, regard- 
less the nature of the interacting region. In calculating 
the current, one needs the spectral function, which is de- 
fined as 

A(tu) = % (G r H - G a H) = i (G>(uj) - G<(uj)) , (7) 

and the lesser Green's function G < , which is proportional 
to spectral function and the occupation of the electron. 

When the operator X is replaced by its expectation 
value, the Hamiltonian Eq.JSJ) is decoupled from the 
phonon operator, then the interacting lesser Green func- 
tion may be separated: 



x=0,T=0 

X=0.6co ,T=0 



>=0.6(D ,T=1 .5cb 



fLiHtq e -iHt 



G< (t) = i (d* (0) d (t)) = % e im d 

\ I el \ I ph 

= G < (t)e-*(- t \ (8) 
and similarly, 

G> (t) = -i (d (i) dt (0)) = G> (t) e-* (t) , (9) 

where G > ^ < - ) (t) is the dressed greater (lesser) Green func- 
tion for a dressed SMT electron governed by H e ; , and the 
factor e - *^*' arises from the trace of the phonon parts 
(Xt(0)X(t)) p/[ or (X(t)Xt(0)) pfc) respectively^ 

$ (t) = g [N ph (1 - e^) + (N ph + 1) (l - e"^ *)] (10) 

Note that because t) ^ &(t), the decoupling approxi- 
mation used in the majority of previous publications ; 16 ! 26 
which directly decouples the retarded (advanced) Green 
function as 



G r[a) {t) = G r{a) (t)e- m , 



(11) 



has ignored the difference between the N p h and N p h + 1 
in the expression of <&(£)• Although such approximation 
works well at high temperature, because N. 



ph 



N, 



ph 



i. 



it does make difference at low temperature because of 
vanishing N p h- We shall further discuss this point latter. 




FIG. 2: The dimensionless Spectral function of the SMT elec- 
tron for different EPI strengths and temperatures. The pa- 
rameters for calculation are taken as: P = 0.2u>o, A*£ = I^R 
• — 0, and the unit Ao = 2/T. The strength of the EPI used 
here is the same as the references-^. 



By using the identity e~*W = T,n=-oc L n e~ mulot , the 
greater and lesser Green functions can be respectively 
expanded as 

oo 

G>H = L n G> (u; ~ ncj ) , 

Ti — — OO 

OO 

G<H = L n G< {uj + ncj ) , (12) 



n— — oo 



where the index n stands for the number of phonons in- 
volved, and L n are the coefficients depending on temper- 
ature and the strength of EPI. At finite temperature, 

L n EE e -9(2N vh + l) eW0 m In (2g^J N ph (N ph + 1)) ,(13) 

where I n {z) is the n-th Bcsscl function of complex argu- 
ment. At zero temperature, L n simply reads 



L n = 



- g ^r n > 



n < 

Thus the spectral function can be expressed as 



(14) 



A(hj) = iL n G > (uj — uuiq) — G < {lu + nwo) '15) 



71 — — OO 



With the help of the equation of motion approach, the 
retarded (advanced) Green function for the dressed elec- 
tron can be analytically evaluated as 



G r(a) (w) = 



1 



lu - e - E r ( Q )(w) 



(16) 
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where the retarded (advanced) self-energy is 



l^fcl 



keL.R 



A(w)=f«T(w). (17) 



For simplicity, in the wide-band limit, both the real and 
the imaginary part of the self-energy, A(w), and T(u>), are 
assumed to be constants independent of energy. Thus, 
when the real part of the self-energy is absorbed into the 
SMT level shift, only the broadening due to the tunnel 
coupling will be considered. With the assumption of sym- 
metric coupling, T L ps = r, the broadening can be 
expressed as f(ui) = (f L +f R )/2 = f . Then the spectral 
function of the dressed SMT electron can be cast into 



A(co) 



2T 



(18) 



Following the Keldysh equation^! for the lesser Green 
function, i.e., G<{uj) = G r (w)£<(u;)G Q (w) with 
E < (cj) = «r(/i(cj) + fn(u)), the relation between the 
dressed greater or lesser Green function and the dressed 
spectral function is respectively as 



G> (u) 



II (w) - f R (u) 



A(u>) 



G (<•>) = * o A H- 



(19) 



Thus the spectral function of the SMT electron A(ui) can 
readily be obtained by substituting Eas. lfH^ and (fl%|) into 
Eq.CES). 



III. THE SPECTRAL FUNCTION OF THE SMT 
ELECTRON 

A. The dressing effects caused by the EPI 

The spectral function for our structure generally de- 
pends on the EPI strength, temperature, and chemical 
potential in two leads relative to the SMT level. Let 
us first address the first two dependencies in the case of 
I^l = = eo- 

As plotted in Fig|2| the effect of the EPI and temper- 
ature on the spectral functions of the SMT electron are 
evaluated. For comparison, the spectral function without 
the EPI at OK is also shown, which exhibits a single res- 
onant peak at eo with a Lorentzian. Compared with the 
EPI-frcc case, with finite EPI, the resonant peak at eo in 
the spectral function is red-shifted by gu>o, and the spec- 
tral peak is sharpened. This results from the renormal- 
ized effects on the SMT level and the dressing effect on 
the tunnelling matrix elements Vk due to the EPI. More 
noticeably, new satellite peaks may appear at Zq ± uloq in 
the spectral function, implying that due to the EPI the 
ground state for the coupled system may possess finite 
components involving n phonons. For later convenience, 
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FIG. 3: The total weights of the +n-th and the — n-th phonon 
sidebands, W n , as functions of the (a) strength of EPI A and 
(b) temperature T, where F = 0.2ljo, while ^l(r) are arbi- 
trary. The dot lines are the numerical summation, W n . 



we label the resonant peak located at eo as the zero- 
phonon-band, and the satellite peaks located at ?o ± n ^o 
as the ±n-th phonon sidebands. In general, these phonon 
sidebands are not symmetric with respect to eo- By using 
the identity 



i I <ko(G > (oj — nojo) — G < (<jJ + nu>o)) 

/>oo 

= i / dcu(G > {oj) -G < {uj)) =2tt, (20) 



and 2~2^ = _ O0 L n = exp{$(0)} = 1, it can be seen that 
the sum rule for the spectral function still holds in the 
presence of EPI, namely 



<kuA(u)) = 2ir. 



(21) 



As shown in Fig'0 the weight of each peak, defined 
as the integrated area under the peak divided by 2ir and 
denoted by W± n , is sensitive to the EPI strength and 
temperature. We can further define sum of the spectral 
weights of the +n-th and —n-th phonon sidebands as 



(e 



na;o/3/2 



-W- r 

-nui j3/2\ 



L, 



-g(2N ph + l) 



)I n [2gJN ph (N ph + 1) (22) 



Fig|3] shows that when A or T increases, the spectral 
weight of the zero-phonon band, Wo decreases mono- 
tonically; while the phonon sideband, W n (n ^ 0), gets 
lager first, then reaches a maximum value for a cer- 
tain EPI strength or temperature, and decreases again. 
This reflects that when increasing the temperature or 
EPI strength, the occupation probability of the phonon 
sideband increases, keeping the conservation of the total 
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FIG. 4: The spectral functions of the SMT electron for four 
typical configurations of Hl(r)'- (a) Hl = Mfl = ?o, (b) 
Mi = eb + 0.9wo, = ?o — 0.5oj , (c) /Lt£ = eg + 2.2cj , fi R = 
eo + 0.8wo, (d) fiL — eo — 0.8cjq, /ir = eb — 2.2wo, where the 
parameters are T = 0.2cjo, A = 0.6tJo, T = and So = 0. The 
balls connected by the dash lines in (a) and (b) denote the 
resonant tunnelling processes for electrons. 



spectral weight, i.e., X^^Lo = 1- Note that these vari- 
ation trends are not in contradiction with Fig|2 where 
the heights of both the zero-phonon peak and phonon 
sidebands increase as temperature gets higher, as the cor- 
responding widths are narrowed simultaneously. 



B. The interplay between the EPI and the SMT 
electron-lead coupling 

Now let us take a close look at the dependence of the 
spectral function of the SMT electron on the chemical 
potentials in two leads in the presence of the EPI. 

As shown in Fig0J it is quite interesting that the spec- 
tral function of the SMT electron profiles quite differently 
for the same EPI strength and temperature (OK) , but un- 
der different configurations of the chemical potentials in 
two leads. And the phonon sidebands, whose lineshapes 
may vary abruptly at certain frequency, can be nonvan- 
ishing in both sides of the zero-phonon band, even at zero 
temperature. This is quite different from the usual Inde- 
pendent Boson Model, according to which, the phonon 
sidebands have nothing to do with the lead chemical po- 
tentials, and at zero temperature, they should vanish on 
the negative energy side of the zero-phonon band of the 
hole. 

The unusual behavior of the present spectral func- 
tion can be understood with Eas. (|15|) and (|19|) . in which 
G > (lo — uloq), and G < (tu + nwo) in general depends on 



the Fermi function of in the left (right) lead, and 

thus on the chemical potential, m.l(M-r)j respectively. At 
zero temperature, the expansion coefficients L n vanish 
for negative n, therefore only the zero-phonon band of the 
spectral function can both the lesser and greater Green 
function contribute to; while the phonon sideband be- 
low the zero-phonon band results from the G < (u> + nu>o), 
and that above the zero-phonon band comes from the 
G > (u> — nwo)- 

Recalling that the lesser and greater Green functions 
G < ' > ) correspond to the dressed SMT electron and hole 
propagator, which are proportional to the occupation 
number, tismt, for the SMT electron, or 1 — ngMT for 
the hole, respectively. Different from the Independent 
Boson Model, where tismt is fixed to either zero or one 
at OK, as the Fermi level for the isolate system is well 
defined, in the present non-equilibrium model, due to the 
tunnelling coupling the occupation of the SMT electron 
is determined by the chemical potentials at the left and 
right leads when the system approaches a steady state. 
Thus, by tuning the lead chemical potentials, Usmt can 
take any value between zero and one, consequently, even 
at OK the contribution from both the electron and hole 
to the spectral function of the SMT need considering. 
At low temperature the SMT electron and hole can only 
emit phonons, so the phonon sidebands below the zero- 
phonon band come from the occupative SMT electron, 
while those above the zero-phonon band originate from 
the SMT hole. If there is a partial occupation in the 
zero-phonon band for the SMT electron, it is required 
to carefully examine the contribution from both electron 
and hole. It is the different treatment of usmt that leads 
to the unusual behavior of the phonon sidebands at zero 
temperature. 

The spectral function can be explicitly expressed as 
functions of the bias and V m id, the average of two lead 
potentials relative to the renormalized SMT level as de- 
fined in the Section II, i.e., A(lo, V m id, Vuas)- In this way 
the effect due to the bias and due to the V m id on the 
spectral function will be considered separately. We have 
found that the V m id mostly affects the lineshape of the 
spectral function in the present model. Four typical con- 
figurations of the lead chemical potentials in Fig0| are 
divided into two categories: (1) The V m id ~ CO) namely 
the SMT level is partially filled. Then for both bias in 
FigEIa) and FigQJb), the resonant tunnelling may take 
place. (2) The V m id deviates from eo significantly com- 
pared with the band-width T, so that the SMT level is 
cither fully occupied or totally empty, i.e. usmt — 1 
(FigQJc)), or tismt — (FigQJd)), in which no resonant 
tunnelling can occur. 

In Fig0{a), both chemical potentials in the leads align 
exactly with the renormalized SMT level, where and the 
lineshape of each phonon sideband exhibits discontinu- 
ity at certain frequency. This is similar to the case 
Ml = fJ-R = e o ( c /- Fig|2J), except for that the spectral 
function in Fig^Ja) is symmetric with respect to eg, while 
that in Figgis asymmetric. This is understandable, be- 
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cause at zero bias and zero temperature there exist a 
well-defined boundary, /i^ = (j, R , separating the SMT 
electron from hole. The symmetric spectral function is an 
exception when the renormalized SMT level happens to 
coincide with the boundary and usmt = 0.5. In FigUfb), 
when the bias is large enough to enclose the most part 
of the zero-phonon band, both distributions of the SMT 
electron or hole in the zero-phonon band become Loren- 
zian, and so does the lineshape of each phonon sideband. 
Strictly speaking, the spectral function in Fig^Jb) is not 
symmetric around 6q as V m id ^ eo> but the asymmetric 
part can be hardly discerned for the present case. It is 
noticed that whenever one of the phonon sidebands en- 
ters the bias region, there will be the phonon-assisted 
tunnelling process. Compare Fig^Jb) to (a), it is ev- 
ident that the bias is not only important in determin- 
ing the profile of each phonon sideband, but also cru- 
cial in controlling the phonon-assisted tunnelling. It is 
interesting to observe that when usmt — (FigEJd)), 
which corresponds to the one-hole picture, the usual re- 
sult following the Independent Boson Model with one 
hole^ is recovered; on the other hand, when usmt — 1 
as in Fig^Jc), which corresponds to the picture of one- 
electron, the spectral function is reversed with respect to 
eo compared with (d). Comparison among the four con- 
figurations shows that the V m id relative to eo has played 
a decisive role in determining the occupation of the SMT 
level, or equivalently the partition between the SMT elec- 
tron and hole. 

The spectral function has been found to have the sym- 
metry as follows, 

A(u, V mi< l, -Vbias) = A(UJ, V mi d, Vbias), (23) 

and 

A(uj = ? - Auj, Vmid = eo/e - AV, V bias ) 
= A(oj = ? + Aw, V mid = ? /e + AV, V bias ). (24) 



Compared to the EPI-free case, the spectral function in 
the EPI is generally asymmetric with respect to ?o unless 
eVnid = £o- This implies a broken electron-hole symme- 
try, which can also be seen from Fig0J where the spectral 
function in (a) is symmetric with respect to eoi while in 
FigOl (b), (c) and (d), there are broken symmetry for 
electron and hole, particularly obvious for cases of (c) 
and (d). 

With increasing temperature, one can expect that the 
profile of the dot spectral function should be less and less 
sensitive to the variation of the lead chemical potentials, 
because the Fermi distribution varies continuously across 
fJ-L(R) a t high temperature, so do the phonon sidebands. 



IV. THE TRANSPORT PROPERTIES 

Based on the spectral function discussed above, in this 
section, we shall investigate the differential conductance 
and tunnelling current of the SMT(QD) system. The case 
of zero temperature will be discussed first, then follows 
the modifications by the finite temperature effect. 



A. The zero temperature 

At zero temperature, the Fermi distribution functions 
are the step functions &{(J-l(r) ~ UJ ), an d the coefficients 
L n reduce to the Eq^J thus the current can be expressed 
explicitly as 



J = Jh ^2 LnT / duj t ~ w ) - ^ R ~ {[6 {hl-oj- nuj Q ) + 9 (p, R - u- nuj )]A(uj + nuj Q ) 

n ** 

+ [2 - 9 (/jll - uj + nuj ) - 9 (fx R - lu + nuo)] A(uj - ra^o)} • (25) 



n=0 



r 



Using ^l(r) = eV m id ± eT4ms/2, and 80(^1^^ — tance can be obtained by performing dJ/dVuas, 
uj)/dVbias = ±ed(fi>L(R) ~ k-0/2, the differential Conduc- 



ed r „ _ _ _ 

G(Vmid, Vbias) = 22 LnT i®( eV bias ~ n^d) A(fJ, L ) + A(/J, R ) + A(fl L - HWq) + A(fl R + HUJ a ) 



n=0 



+@(-eV bias - nu Q ) A(fj, L ) + A(fj, R ) + A(fj, L + uujq) +A((j, r - uujq) i. 



(26) 
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FIG. 5: (Color on line) The dimensionless differential con- 
ductance (Go = e 2 /2h) at T = OK as functions of V m id and 
Vbias- The parameters for calculation are taken as: T = 0.2ljo, 
eo = wo, and eo = which implies that for the sake of clarity 
a stronger EPI strength, A = Uo, is chosen. The color scale 
runs from zero (blue) to Go (red). The above (a) and right 
(b) panels to the map are the sections for a fixed value of 
e ^Mas = 3-6^o and eV^ii = 0.8ljo, respectively, marked on 
the map with dash lines. 

It is easy to verify that the tunnelling current and differ- 
ential conductance satisfy the symmetry relations as 

J(V mid = e /e - AV, V b%as ) = J(V mid = e /e + AV, V Uas ), 
G(V mid = e /e - AV, Vitas) = G(V mtd = e /e + AV, V bias ), 

and 

J{Vmid, Vbias) J (Vmid , Vbias ) , 

G (Vmid,— Vbias) = G (V m id, Vbias)- (27) 

It is noticed that due to the equivalence between the 
current of electron type and hole type, the broken sym- 
metry with respect to eo in the spectral function of SMT 
electron is now restored for the tunnelling current and 
differential conductance. 

The map plotted in Figgis the calculated differential 
conductance as function of V m id and Vuas in the pres- 
ence of the EPI. The dashed lines in FigOJa) and (b) 
represent the differential conductance for a fixed value of 
eViias = 3.6wo and eVj™ d = 0.8w , respectively; while 
the solid lines correspond to the differential conductance 
without the EPI as for reference, in which the peaks ap- 
pear only when one of the lead chemical potentials is 
aligned with the SMT level, i.e., (j,l((J-r) = eo- Com- 
pared with the EPI-free case, the corresponding peaks 
in the presence of finite EPI are sharpened, and red- 
shifted by gtuo, and more noticeably, new satellite peaks 
are shown up, which are associated with the phonon- 
assistcd tunnelling processes as depicted in FigQJb), 



FIG. 6: (Color on line) The map of tunnelling current, 
J (Vmid, Vbias), at T = OK, where the parameters are taken 
as the same as the Fig|KJ and Jo = eF /2H. The color scale 
runs from —0.9 Jo (blue) to 0.9 Jo (red). The solid and dashed 
lines in sections (a) and (b) represent the tunnelling current 
for the cases without and with EPI, respectively. 



when (i L = e + nuj or fj, R = e - nu Q for V blas > 0; 
Mz, = eo - nuj or fx R = e + nuj for V bia s < 0, where 
n < Q(\eV bias \/u!o). 

The map plotted in Figgis the tunnelling current as 
function of Vmid and Vbias, which can be divided into 
several plateaus with the bounded lines corresponding to 
the differential conductance peaks in Fig|SJ It should be 
pointed out that the phonon-assistcd tunnelling can take 
place even at zero temperature. Although no thermal 
phonon is available at OK, the phonon-emitted process 
accompanying the tunnel is possible, if the phonon en- 
ergy can be supplied by the bias voltage across the SMT. 
As also shown in the map (FigJHJ, the phonon-assisted 
peaks are absent in the left and right quarters, which are 
bounded by the lines of the zero-phonon peaks. 

Now let us discuss the relationship between the trans- 
port properties and the SMT spectral function. As men- 
tioned above, due to the dressing effects of EPI, the SMT 
electron still has a finite probability to occupy the phonon 
sidebands even at zero temperature. Once a phonon side- 
band of the SMT enters the bias region, the phonon- 
assistant-channel would be open, as depicted in FigQJb). 
Thus the information of the SMT spectral function can 
be inferred from the spectra of the tunnelling current or 
differential conductance. For example, for a fixed V m id, 
the integral of the differential conductance over Vbias sat- 
isfies a sum rule, dVbiasGiVmid, Vuas) = Jo, which 
just results from the conservation of the spectral weights 
in Ea. l|21|l . While for a fixed Vbias, the integral of the 
differential conductance over V m id does not equal Jo in 
consequence of the transfer of a finite spectral weights 
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C./G, 




5 -1.5 -3 



FIG. 7: (Color on line) The dimensionless differential con- 
ductance and tunnelling current as functions of V m id and 
V bias , (a) G(V mid , Vuas), (b) J{Vm id , V b ias), where V = 0.2cj , 
A = ljo, T = 0.6tJo, eo = 0, and the units Go = e 2 /2h, 
Jo = eT/2h. The color scale runs from zero conductance 
(blue) to 0.05Go (red) in (a) and from —0.15 Jo (blue) to 
0.15 Jo (red) in (b). 



Green function is decoupled according to Ea. (|llf) . the 
SMT spectral function would be expanded in the follow- 
ing way, 



B(u) 



E 



(28) 



At zero temperature, we have n > 0, the profile of B(u>) is 
similar to Figged), namely there is a set of non- vanishing 
phonon sidebands above the zero-phonon band, which 
is independent of the Fermi distributions in the leads. 
Thus, following the approximation i|ll|) . in the zero bias 
limit, the differential conductance would exhibit a set of 
phonon peaks. Physically, our results which follow Eqs. 
@ and (JSJ, seem more reasonable, because the tunnelling 
electron can neither absorb nor emit any phonon in the 
zero temperature and zero bias limit. 



to those phonon sidebands, which do not participate in 
tunnelling. This is a little different from the tunnelling 
process without the EPI, where the sum rule for the dif- 
ferential conductance always holds. 

We will end this subsection with a comparison between 
the results based on different decoupling approximations, 
i.e., Eq.© and Eq. Hll|) . In the zero bias limit, by us- 
ing Eq.©, the differential conductance would have no 
phonon-assisted peaks as clearly shown in Fig|Hj On the 
other hand, by following Ea. (|ll|l . there would be a set 
of non-vanishing phonon-assisted peaks in the differen- 
tial conductancei£*i§i2£i2£i2S. This discrepancy can be 
explained in terms of the different SMT spectral func- 
tions obtained in different approximations. Once the 



B. finite temperature 

The tunnelling current at finite temperature can read- 
ily be obtained from Eq.©. With the help of the 
identities L_„ = e~^ nu3 °L n and d f L ( R) {uj) / dV bi as = 
±e/3/i(fl)(w)(l — /i(H)(w))/2, the differential conduc- 
tance can be expressed in a compact form: 

g2 p oo -oo 

G=— 22 L n duF n {w)A{u - nu ), (29) 

n— — oo 00 

where the factors F n (u>) depend on the Fermi distribu- 
tions in the leads through, 



F n (u>) = ? \ [/lH(1 - /7») + /k(w)(1 - /aM)][1 + 2 - ™ ) + f R {u> - nw ))] 

- fn(u)))[f L (u; - nu ){l - - nuo)) - /k(w - nw )(l - /r(w - nui ))} 

I 



2 

e -/3no) _ I 



Both maps for the differential conductance and tun- 
nelling current at finite temperature are plotted in Fig0 
Compared to the zero temperature case, the differential 
conductance peaks at finite temperature are broadened 
and smeared out to some extent, and the tunnelling cur- 
rent profile become smoother and less sensitive to the 
bias, which results from the smoother Fermi distributions 
in the leads at finite temperature. Moreover, since the 
phonon sidebands of the spectral function will distribute 
more symmetrically around the zero-phonon-peak at fi- 
nite temperature than the OK case as shown in FigEl 
the phonon-assisted peaks of the differential conductance 
might appear in the two forbidden quarters as shown in 
FigEJat OK. Therefore, the difference caused by the dif- 



ferent decoupling approximations as discussed above will 
be diminished with increasing the temperatures. 



V. CONCLUSION 

In summary, we have theoretically studied the non- 
equilibrium quantum transport through the single 
molecule transistor or quantum dot in the presence of 
the local electron-phonon interaction. Owing to the EPI, 
the spectral function of the dot electron will manifest it- 
self in the phonon-dressing effects, such as the red-shift 
and sharpening of the zero phonon peak, and the transfer 
of finite spectral weights to the phonon sidebands. Fur- 
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thcrmore, due to the interplay between the EPI and the 
dot-lead coupling, the zero-phonon band and the phonon 
sidebands of the spectral function at low temperature 
critically depend on the chemical potentials at two leads 
and can thus be manipulated. Although the sum of the 
spectral weights for the +n-th and — n-th phonon side- 
band is fixed for a given temperature and EPI strength, 
the distribution of the spectral weights is dependent on 
the chemical potentials in two leads, in particular on the 
average of two lead potentials relative to the renormal- 
ized SMT level, showing the broken electron-hole symme- 
try. The tunnelling current and differential conductance 
have been analyzed at both zero and finite temperatures, 
which also reveal the dependence on the chemical po- 
tentials of the leads. Different approximations used in 



decoupling the electron-phonon system have been com- 
pared and discussed. It has been found that although 
different decoupling approximations can yield nearly the 
same results at high temperature, they do lead to quite 
different behaviors at zero temperature. 
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